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C/^ ■ Abstract Based on the turbulent convection model (TCM) of Li & Yang (2007), we have 

(— I I studied the characteristics of turbulent convection in the envelopes of 2 and 5Mq stars 

Qh' at the RGB and AGB phases. The TCM has been applied successfully in the whole con- 

JL ' vective envelopes including the convective unstable zone and the overshooting regions. 

^ . We find that the convective motions become stronger and stronger when the stellar mod- 

'^ • els are located upper and upper along the Hayashi line. In the convective unstable zone, 

Ci ', we find that the turbulent correlations are proportional to functions of a common factor 

(V — y ad)T, which explains similar distributions of those correlations. For the TCM we 

find that if the obtained stellar structure of temperature is close to that of the MLT, the 

^ . convective motion will be at a much larger velocity and thus more violent. However, if 

f^ ' the turbulent velocity is adjusted to close to that of the MLT, the superadiabatic convection 

\^ , zone is much more extended inward and a lower effective temperature of the stellar model 

C^^ ' will be obtained. For the overshooting distance we find that the e-folding lengths of the 

turbulent kinetic energy k in both the top and bottom overshooting regions decrease as the 

f~~. ■ stellar model is located up along the Hayashi line, but both the extents of the decrease are 

^^ , not obvious. And the overshooting distances of the turbulent correlation u'^T' are almost 

the same for the different stellar models with a same set of TCM's parameters. For the 

decay way of the kinetic energy k we find that they are very similar for the different stellar 

models based on a same set of TCM's parameters, and there is a nearly linear relationship 

between Ig k and In P for the different stellar models. When Cg or a increases while the 

other parameters are fixed, the obtained linearly decaying distance will become longer. 
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1 INTRODUCTION 

Red giant branch (RGB) and asymptotic giant branch (AGB) stars are characterized by the convective 
motion in their envelopes, which can extend over an enormous range in mass (or radius). Because of 
huge scale of convection zones and small viscosity of the stellar material, turbulent flow instead of 
laminar one always occurs in their outer envelopes. Turbulent convection may result in many important 
effects in stars: mixing different elements to be homogeneous in the convection zone and adding fresh 
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fuel to the nuclear burning region, dredging the internal material processed by H or He burning up to the 
stellar surface, and being as an important means for heat transfer These effects significantly influence 
the structure and evolution of stars. Many observational phenomena appearing in the RGB and AGB 
stars are ascribed to the convection. For the RGB stars as an example, many of them show chemical 
anomalies (Pilachowski et al. 1993; Charbonnel 1995; Gratton et al. 2000; Kraft et al. 1993; Ramiez & 
Cohen 2002; Gratton et al. 2004; Busso et al. 2007; Recio-Blanco & de Laverny 2007). And we are still 
perplexed to the so-called carbon star mystery (Iben 1981; Iben 1975, 1977; Sackmann 1980; Lattanzio 
1989; Hollowell & Iben 1988; Straniero et al. 1997; Herwig et al. 1997; Herwig 2005) for the AGB stars. 
Furthermore, the observed mass loss occurring in the RGB and AGB stars may be due to the turbulent 
pressure (Jiang & Huang 1997). To solve these problems we need to accurately deal with the convection 
and to learn details of its characteristics. 

For some RGB and AGB stars the convective motion near the surface of the stellar envelope are 
superadiabtic because of less efficient convective energy transport, and may sometimes become super- 
sonic based on the mixing-length theory (MLT) (Deng & Xiong 2001). The MLT (Bohm-Vitense 1953, 
1958) is commonly used to deal with the convection. However, due to its simplicity many shortcomings 
are found (Renzini 1987; Baker & Kuhfuss 1987; Spruit et al. 1990; Pederson et al. 1990), which thus 
lead to some confused and debatable problems. For example, convective overshooting (e.g. Saslaw & 
Schwarzschild 1965), semiconvection (e.g. Castellani et al. 1971; Sreenivasan & Wilson 1978) and the 
breathing convection (e.g. Castellani et al. 1985). Therefore new convection models have been proposed 
successively, for example, non-local mixing length theories (Ulrich 1970; Maeder 1975; Bressan et al. 
1981) and full-spectrum convection theory (Canuto & Mazzitelli 1991). However, they are still based 
on the framework of the MLT. A better theory to overcome the MLT's drawbacks is the turbulent con- 
vection models (TCMs), which are deduced directly from the Navier-Stokes equations, and can describe 
many turbulence properties (Xiong 1979, 1981, 1985, 1990; Canuto 1992, 1994, 1997; Xiong, Cheng 
& Deng 1997; Canuto & Dubovikov 1998; Canuto et al. 1996). However, the dynamic equations for 
turbulent correlations are expressed by higher order correspondents due to the inherent non-linearity of 
the Navier-Stokes equations. Therefore these equations are needed to be cut-off with some closure ap- 
proximations in order to be applicable in calculations of the stellar structure and evolution. Nevertheless, 
some free parameters have to be introduced by the closure assumptions. Different closure assumptions 
correspond to different TCMs, and many of them may not be easy to be incorporated into a stellar evolu- 
tion code. Recently a simple TCM proposed by Li & Yang (2007) was successfully applied to the solar 
models (Yang & Li 2007), in which the introduced parameters of the TCM are directly related to the 
corresponding turbulence physics. Some significant changes in the structure of the solar convection zone 
and better results compared to the solar p-mode observations are obtained. Shortly afterwards, Zhang & 
Li (2009) applied it successfully to the overshooting region of the solar convection zone. 

Up to now, the TCMs have seldom been applied to very large convective envelopes of stars, e.g. in 
the RGB or AGB stars, which are rather different from the solar environment. Furthermore, in view of 
the key roles of convection played in these stars, we choose three sites of stellar evolution to test the 
TCM of Li & Yang (2007). We try to find out convection characteristics for two stars of 5 and 2Mq 
at the RGB and AGB phases, and to find out their dependence on the TCM's parameter. To analyze 
the properties of the turbulent correlations among velocity and temperature fluctuations, we separate 
the whole convection envelope (derived self-consistently from the TCM) into the convective unstable 
zone and the convective overshooting regions, respectively. We firstly describe the basic equations of the 
TCM in Section 2. The information of the evolutionary code and input physics are described in Section 
3. In Section 4, the convection characteristics in the convective unstable zone and overshooting regions 
are presented and discussed, respectively. The dependence of the structure of the overshooting regions 
on the TCM's parameters is given in Section 5. Finally, some concluding remarks are summarized in 
Section 6. 
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2 EQUATIONS OF TCM 

Using the Reynolds decomposition approximation, and with the aid of some closure approximations, 
the second-order moment equations of the TCM are (Li & Yang 2007 for details); 
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The terms on the left-hand side of Eqs.(l)-(4) are diffusion ones of the turbulence, representing 
the non-local characteristics of turbulence and describing mainly the properties of the overshooting 
regions outside the Schwarzschild boundaries of a convection zone. The terms on the right-hand side 
of Eqs.(l)-(4) are some production and dissipation ones of the turbulence. All the symbols appeared in 
above equations have their usual meaning (Li & Yang 2001; Li & Yang 2007). There are seven free 
parameters in above equations. They are Ct, Ce, Ck, Cti, Cei, Cg, and a. The former three ones are 
dissipation parameters which describe respectively the dissipation of the turbulent correlation u^T', the 
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dissipation of the temperature fluctuation T , and actually the anisotropic degree of the turbulence. The 
larger the values of them are, the smaller the convective heat fluxes will be (Li & Yang 2007). The next 

three parameters are diffusion ones related respectively to the three turbulent correlations u^T', T' 
and u'^u'^. Larger values of them will result in lowered and expanded profiles of these correlations (Li 
& Yang 2007). The last parameter a, the ratio of a typical length I to the local pressure scaleheight, is 
introduced by the closure assumption similar to the mixing length parameter in the MLT. The effects of 
these parameters on the properties of turbulence, and then on the structure of the solar convection zone, 
are investigated by Li & Yang (2007) and Zhang & Li (2009). Here we choose several sets of these 
parameters' values to study the properties of the turbulent convection in the RGB and AGB stars. 



3 INPUT PHYSICS 

We calculate the evolution of two stars of 2 and SM© respectively. The initial composition (Z=0.02, 
X=0.7) is fixed for both stellar models. 

The free parameters adopted in the TCM can be determined by the fluid dynamics experiments in 
the terrestrial environment (e.g. Hossain & Rodi 1982). Given a set of appropriate values of dissipation 
and diffusion parameters, the value of a can be determined by calibrating the solar model (Li & Yang 
2007; Zhang & Li 2009). The value of a obtained in this way can be used to RGB or AGB calculations, 
but this approach is at the sacrifice of the uncertainty of its value (Herwig 2005). In this paper, we adopt 
several sets of the parameters' values, which are derived directly from or mainly based on the works 
of Li & Yang (2007) and Zhang & Li (2009) for the solar model and are shown in Table 1 . In which 
the values of Ct, Ce and Ck are mainly derived from the terrestrial experiments and their typical values 
are 3.0, 1.25 and 2.5, respectively. Nevertheless, Ct = 7.0 is taken from the work of Canuto (1993) 
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Table 1 The information of the TCM's parameters 



Model Name 


Ct 


a 


Ck 


Cti 


Cel 


c. 


a 


MLT 














1.70 


LM 


3.0 


1.25 


2.5 








0.90 


NLMa 


3.0 


1.25 


2.5 


0.03 


0.03 


0.03 


0.90 


NLMbl 


7.0 


0.20 


2.5 


1.0 X 10^^ 


1.0 X lO"'^ 


0.10 


0.76 


NLMb2 


7.0 


0.20 


2.5 


1.0 X 10^^ 


1.0 X lO"'^ 


0.10 


0.74 


NLMb3 


7.0 


0.20 


2.5 


1.0 X 10"^ 


1.0 X 10"^ 


0.10 


0.15 


NLMcl 


3.0 


1.25 


2.5 


0.05 


0.05 


0.05 


0.68 


NLMc2 


3.0 


1.25 


2.5 


0.05 


0.05 


0.04 


0.68 


NLMc3 


3.0 


1.25 


2.5 


0.05 


0.05 


0.03 


0.68 


NLMc4 


3.0 


1.25 


2.5 


0.05 


0.05 


0.02 


0.68 


NLMc5 


3.0 


1.25 


2.5 


0.05 


0.05 


0.01 


0.68 


NLMd 


0.2 


0.10 


2.5 


0.15 


0.25 


0.10 


0.05 



and Ct = 0.20, Cg = 0.10 from the suggestion that smaller values of them will lead to a better result 



compared to the solar p-mode observation (Yang & Li 2007). Furthermore, Cti = Cei 
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adopted, meaning a negligible diffusion effect of u'^T' and T'^ and a local convection model of them, 
to assure that an almost same temperature structure in stellar model can be obtained for both the MLT 
and TCM and then a meaningful comparison between them can be done, which will be discussed in 
the following subsection 4.3. For simplicity we call the TCM with the diffusion terms in Eqs.(l)-(4) 
as NLMs and the TCM without the diffusion terms as LM. The MLT with a = 1.70 is also adopted 
to calculate the stellar convection at the RGB and AGB phases for comparisons. In Table 1 both the 
parameter sets of NLMa and NLMd are corresponding to the standard solar model calibration. 

The stellar evolution code used in the present paper was originally described by B. Paczynski and 
M. Kozlowski and updated by R. Sienliewicz. Nuclear reaction rates are adopted from BP95 (Bahcall, 
Pinsonneault & Wasserburg 1995). The equation of state is the OPAL EOS from Rogers (1994) and 
Rogers, Swenson & Iglesias (1996). The OPAL opacities GN93hz series (Rogers & Iglesias 1995; 
Iglesias & Rogers 1996) are used in the high-temperature region. In the stellar atmosphere, low- 
temperature opacities from Alexander & Ferguson (1994) are used. 

To study the convection properties based on the TCM, we select three evolved models during the 
RGB and AGB phases derived from the MLT, which are indicated by solid stars in Fig.l. 




Fig. 1 The HR diagram for the stellar models of 2 and 5Mq based on the MLT. The 
three solid stars indicate the certain locations of the stellar models during the RGB 

(2Mq : IgL/Lg = 1.58, Te// = 4750K;5Mq : IgL/^o = 3.09, Te// = 4190/i) 
and AGB (SAf© : Igi/Lo = 3.81, T^jf = 3570A') phases 
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Fig. 2 Distributions of the turbulent corre- 
lations for the 5A/0 star at the AGB phase. 
The solid lines are for NLMa, dashed lines 
for LM. 



Fig. 3 The same as Fig. 2, but for the 5AIq 
star at the RGB phase. 



4 CHARACTERISTICS OF TURBULENCE 



4.1 Profiles of tlie Turbulent Correlations 



The profiles of the turbulent correlations vfc, u'^T' and T'^ are displayed in Fig. 2 for the 5Mq model at 
the AGB phase. It can be seen that convection widely develops in the stellar envelope (3.5 < Ig T < 6.4) 
and the turbulent kinetic energy k maintains at a high level in most of the convection zone. Particularly, 
there are several peaks for all the three correlations appearing at almost the same locations. The ap- 
pearance of the maximum turbulent kinetic energy is mainly due to ionization of hydrogen and helium 
which intensively blocks the transfer of heat and results in strong buoyancy to drive the convective mo- 
tion. The profiles of the correlations u'^T' and T'^ vary in a similar way as that of Vfc, except that the 
largest peaks are much higher than that of vfc. 

Diffusion terms in the turbulence equations are taken into account in NLMa but not in LM. As a 
result the non-local effect is behaved remarkably for NLMa around the boundaries of the convective 
unstable zone, especially around the upper boundary where the NLMa's solution significantly deviates 
from the LM's. However, the non-local effect can be safely ignored in the interior of the convective un- 
stable zone because the turbulent correlations distribute quite homogeneously there. These conclusions 
agree well with those of Yang & Li (2007) and Deng & Xiong (2008). 

However, for the other two stellar models of 5 and 2Mq at the RGB phase, the profiles of the three 
correlations, which are respectively shown in Fig. 3 and Fig.4, have some differences for both NLMa 
and LM. These two models have smaller values of the three correlations. It means that the turbulent 
motion is less violent and thus transfers less heat. Furthermore, except for the maximum peak of those 
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Fig. 4 The same as Fig. 2, but for the 2Mq star at the RGB 
phase. 

correlations the other peaks become less prominent, especially for the model of 2Mq at the RGB phase. 
Besides, near the bottom of the convection envelope the turbulent kinetic energy k has a weak peak, and 
its value seems to decrease monotonically with the stellar luminosity. In another word, the convective 
motion becomes stronger and stronger when the stellar models are located up along the Hayashi line. 

4.2 Turbulence Properties in the Convective Unstable Zone 

In order to understand the properties of turbulence in the convective unstable zone, some assumptions 

are invoked to obtain the explicit relationships between the correlations Vk, u[.T' and T' . According 
to the result stated in the former section, the diffusion terms on the left hand sides of Eqs.(l)-(4) can be 
approximately ignored in the convective unstable zone. With the aid of Eqs.(l)-(4) we get: 



rp/2 __ 



where D ~ iS^ 



Eq.(5) can be further reduced to be; 
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By utilizing Eq.(6) and Eq.(2), the velocity-temperature correlation can be approximated as: 



(5) 



(6) 



(7) 



(8) 



It can be seen that all the three correlations depend on (V — Vad)T, which explains why the three 
correlations peak at almost the same places in the convective unstable zone. Fig. 5 shows the profiles 
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of k,KTf' '~'^"" 



[T'^) and (V - ^ ad)T for NLMa of the SM© star at the AGB phase. It can be 



2/3 



1/2 



seen that the profiles of {u'^T') ' , (T'^) and (V — V ad)T behave in a very similar way in most 
of the convective unstable zone. However, the profile of k deviates significantly from the others, but its 
overall trend is almost consistent with them. It should be noted that the above results are suitable for 
other stellar models. 
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Fig.5 Profiles of fc, {u'XY'", (T'^) 
5M0 star based on NLMa. 



1/2 



and (V - Vad)T at the AGB phase of the 



4.3 Comparisons with the MLT 
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Fig. 6 The temperature gradient (including the real, radiative and adiabatic ones) and the 
square root of the turbulent energy vfc as functions of the temperature for the SM© star at the 
AGB phase (the same as in Fig. 2). The dashed line is for MLT, the solid lines in left panels 
are for NLMbl and the solid lines in right panels are for NLMb3. 
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Results of two TCMs (NLMbl and NLMb3) are compared with that of the MLT. Firstly, for the 
comparison between NLMbl and the MLT in the left panel of Fig. 6, we can find that the temperature 
gradient V are almost the same with each other, which means that NLMbl will result in a similar stellar 
structure as the MLT does in the convective unstable zone. However, the value of the turbulence velocity 
(^ \/k) of NLMb 1 is much larger than that of the MLT. In another word, for a fixed heat flux transferred 
by the convective motion, a larger value of convective velocity is requested for the TCM. This is partly 
because that the dissipation effects of the turbulent convection, measured by the parameters Ct, Ce and 
Cfe, are fully taken into account. 

When decreasing the value of parameter a, for example a = 0.15, however, we obtain an almost 
equivalent profile of the convective velocity as the result of the MLT, which is shown in the right panel 
of Fig. 6 (namely NLMb3). It can be found that the temperature gradient V of NLMb3 is much larger 
than that of the MLT near the top of the convection envelope and extending to IgT « 4.7 where it 
approaches the adiabatic temperature gradient Wad, a place much deeper than the case of the MLT 
being approximately at IgT w 4.1. It means that the superadiabatic convection region is much more 
extended inward for the TCM than the case of the MLT and the stellar models based on the TCM will 
be of lower effective temperature than that of the MLT. 

The results obtained above are mainly based on the star of 5Mq at the AGB phase. Likewise, these 
conclusions are still valid for the 2Mq star at the RGB phase and the 5Mq star at the RGB phase. Fig. 7, 
as an example, shows the case of the 2Mq star at the RGB phase. It should be noted from the right panel 
of Fig. 7 that the superadiabatic convection region does not extend so much inward as in the case of the 
BMq star at the AGB phase. 
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Fig. 7 The same as Fig. 6, but for the 2Mq star at the RGB phase. 



4.4 Turbulence Properties in the Overshooting Region 



The turbulent correlations in the top and bottom overshooting regions are shown in Figs. 8-10 for the 
considered stellar models. It can be found that the turbulent kinetic energy k monotonously decreases 
outwards and tends to be zero from the boundaries of the convective unstable zone into the overshooting 

regions. For the correlation T'", there is a peak in both the top and bottom overshooting regions, its 
value being positive and much larger than that of Zhang & Li (2009). There is an exception in the top 
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overshooting region for the 2A/q star at the RGB phase seen in the left panel of Fig. 10, the peak of T'^ 
disappearing completely. Nevertheless, the value of the correlation u'^T' is always negative in both the 
top and bottom overshooting regions. As a result, the temperature gradient V is greater than the radiative 
temperature gradient Vr but smaller than the adiabatic one Vad- 
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Fig. 8 Distributions of s/k, T'^ and the tem- 
perature gradients in the top overshooting re- 
gion (TOV) and the bottom overshooting re- 
gion (BOV) of the 5A/0 star at the AGB 
phase according to NLMa. 



Fig. 9 The same as Fig. 8, but for the 5Mq 
star at the RGB phase. 



On the other hand, the overshooting extents are different for the different correlations in a certain 
overshooting region. The profile of \fk approximately demonstrates the distance that a convective cell 
can reach beyond the boundaries of the convective unstable zone. Its overshooting distance and decaying 
way will be very important to the chemical mixing in the overshooting region, if it is used to construct 
the diffusion process of chemical mixing (e.g. Freytag et al. 1996, Herwig et al. 1997, Salasnich et al. 
1999, Ventura & D'Antona 2005). For the three stellar models, our numerical results show that the e- 
folding length of \fk in the top overshooting region decreases as the stellar model is located upper and 
upper along the Hayashi line (they are 0.37, 0.34 and 0.34, respectively, in the unit of the local pressure 
scaleheight Up), and the case in the bottom overshooting region behaves in the same way (they are 0.24, 
0.23, 0.21i7p, respectively), but the difference among them is not obvious in both the top and bottom 
overshooting regions. However, for the case of u^T' that can be inferred by the length of the region 
satisfying Vad > V > Vi?, all the overshooting distances in the bottom overshooting regions for the 
three models are about 0.15i?p. It means that the bottom convective overshooting has a similar effect 
on the stellar structure for all of the stellar models. In the top overshooting regions, however, it can be 
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seen that V and Wr nearly completely overlap with each other, except for a tiny difference for the 2Mq 
star at the RGB phase. This may be due to very low density in the top overshooting regions, resulting in 
the convective heat flux comparatively rather low. 
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Fig. 10 The same as Fig. 8, but for the 2A/q star at the RGB phase. 



5 DEPENDENCE OF THE STRUCTURE OF THE OVERSHOOTING REGIONS ON TCM'S 
PARAMETERS 

The dependence of convection characteristics of the sun on the TCM's parameters has been investigated 
extensively (Li & Yang 2007; Zhang & Li 2009). The turbulent correlations are mainly affected by the 
corresponding dissipation and diffusion parameters. The ratio of the radial kinetic energy to the hori- 
zontal one u'j.u'^ / u'f^u'i^ are shown in Fig. 11 for the AGB model of 5AIq, focusing on the effect of value 
of Cs on the convective structure in the overshooting regions, because Cg is the most sensitive param- 
eter to determine the overshooting distance (Zhang & Li 2009). It can be seen that with the convective 
cells penetrating more and more into the top overshooting region the convection becomes more and 
more horizontally dominated. Particularly, the turbulence is isotropic when u'^u'^, / u'/^u'i^ = 0.5. It can be 
seen that the larger the value of Cs is, the longer the isotropic turbulence will be developed in the top 
overshooting. With respect to the bottom overshooting region, the ratio u^m^/m^m^ shows a much more 
complicated behavior for smaller values of Cs- When Cs is large enough, this complication disappears 
and the convection becomes more and more horizontally dominated as the convection penetrates inward 
into the overshooting region. 

The decaying law of the turbulent kinetic energy fc is a crucial factor determining the efficiency of 
convective mixing in the overshooting regions. The decay way of the turbulent kinetic energy k in the 
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Fig. 11 Ratio of the radial turbulent kinetic energy to the horizontal one u^,u^,/u^u^j with 
different TCM's parameters for the SA/q star at the AGB phase. 
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Fig. 12 The decay way of the turbulent energy k in the top overshooting zone (TOV) 
and the bottom overshooting zone (BOV) of the 5Mq star at the AGB phase according 
to some non-local models. 
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Fig. 13 The same as Fig. 12, but for the 5Mq star at the RGB phase. 



overshooting regions are given in Figs. 12-14 for the three stellar models based on different choices of 
the TCM's parameters. There is a nearly linear relation between Ig k and In P before a turning point on 
the curve of Ig k, beyond which k begins to drop promptly to zero. For example, for NLMa the relation 
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Fig. 14 The same as Fig. 12, but for the 2Mq star at the RGB phase. 



can be approximated as Ig fc = 2.721nP + Cinthe top overshooting regions, where C = —16.8,-12.4 
and —10.2, respectively, for the three stellar models up along the Hayashi line. On the other hand, in 
the bottom overshooting regions we can obtain similar relations, for instance, for NLMc3 the result 
is Igfc = -4.98 InP + C, where C = 161, 133.8 and 140 for the three stellar models. It should be 
noted that the slope of the decaying law seems only a function of the TCM's parameters, and the stellar 
parameters only affect its constant. 

The extension of these linearly decaying region and the slope of the decaying law are sensitive to 
the TCM's parameters. It can be found that when Cg increases while the other parameters are fixed (e.g. 
NLMcl, NLMc3, NLMc5), the turbulent kinetic energy decays slower and slower, and the obtained 
linearly decaying distance will become longer. For example, for the model of 5Mq star at the AGB 
phase in Fig. 1 2, the length of linearly decaying distance in the bottom overshooting region for NLMc 1 , 
NLMc3, and NLMc5 are respectively about 1.2Hp, lAHp and 1.7Hp. On the other hand, when the 
value of parameter a increases, the slope of the decaying law becomes less and the linearly decaying 
distance is significantly prolonged as seen by comparing the result of NLMd and that of NLMa. From 
Figs. 12-14 it can be found that the decaying ways of the k are very similar for stars of different masses 
and evolutionary stages with same values of the TCM's parameters. 



6 CONCLUDING REMARKS 

Based on the TCM of Li & Yang (2007) we have obtained the characteristics of the turbulent convection 
in the RGB and AGB stars with huge convection zones in their envelopes. Some sets of suitable values 
of the TCM's parameters that have been used in the solar convection zone are adopted to make an 
exploratory application in two stars of 2 and 5Mq at the RGB and AGB phases. In practical calculations 
the TCM has been applied to the whole convective envelope including the overshooting regions. When 
analyzing the characteristics of the turbulent correlations, we separate the whole convective envelope 
into the convective unstable zone and the convective overshooting region. In the convective unstable 
zone, we find an approximation approach to associate the turbulent correlations explicitly, and can be 
used to explain the turbulent properties. Later, we make a comparison between the results of the TCM 
and those of the MLT on the velocity profile of turbulence and in influence on stellar structure and 
evolution. These obtained features of convective motion by the TCM are highly sensitive to the TCM's 
parameters. The main conclusions are summarized as follows: 

1) The non-local effect of the turbulent convection is only exhibited significantly around the bound- 
aries of the convective unstable zone, especially around the upper boundary, which is in agree- 
ment with the results of Li & Yang (2007). The three correlations y/k, u',.T' , and T'^ are nearly 
peaked at the same place in the convective unstable zone for a same set of the TCM's parameters. 
It is because that all of them are directly correlated to a function of (V — Vq^JT. 
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2) For a fixed heat flux (closely related by u'^T') transferred by the convective motion, compared to 
the results of the MLT, the TCM will result in a larger turbulent velocity and thus more violent 
convective motion. Furthermore, for a nearly same profile of the convection velocity obtained by 
adjusting the value of the mixing length parameter a for both the TCM and MLT, the superadia- 
batic convection zone will be extended much more inward for the TCM than for the case of the 
MLT and the stellar models based on the TCM will be of lower effective temperature. 

3) For the stellar models located up along the Hayashi line, we find that the convective motions 
become stronger and stronger. However, the e-folding lengths of vfc in both the top and bottom 
overshooting regions decrease as the stellar model is located up along the Hayashi line. The 
overshooting length of u'^.T' is found to be shorter than other turbulent quantities for a same set 
of the TCM's parameters in both overshooting regions, which is consistent with the results of 
Deng & Xiong (2008) and Zhang & Li (2009). And the overshooting distances of m^T' in the 
bottom overshooting region are almost the same for different stellar models with a same set of 
the TCM's parameters. Therefore, the bottom convective overshooting has a similar effect on the 
stellar structure for all of the considered stellar models. 

4) The diffusive parameter Cg plays an important role in the development of the isotropic turbulence 
in the top overshooting region, and the longer the isotropic turbulence will be obtained for larger 
Cs- In the bottom overshooting region, the convection becomes more and more horizontally 
dominated for larger value of Cs, and the convection shows a much more complicated behavior 
when Cs is small enough. 

5) The decaying ways of the turbulent kinetic energy k are very similar for the different stellar 
models based on a same set of TCM's parameters. We find that there is a nearly linear relation 
between Ig k and In P in most of the overshooting regions. The slope of the decaying law seems 
only a function of the TCM's parameters. The values of Cs and a have a similar effect on the 
decaying distance of the turbulent kinetic energy, and the larger the values of them the slope of 
the decaying law becomes smaller and therefore the linearly decaying distance becomes longer. 
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